Halo Geometry and Dark Matter Annihilation Signal 
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We study the impact of the halo shape and geometry on the expected weakly interacting massive 
particle (WIMP) dark matter annihilation signal from the galactic center. As the halo profile in 
the innermost region is still poorly constrained, we focus on geometrical distortions and consider 
different density behaviors like flat cores, cusps and spikes. We show that asphericity has a strong 
impact on the annihilation signal when the halo profile near the galactic center is flat, but becomes 
gradually less significant for cuspy profiles, and negligible in the presence of a central spike. However, 
the astrophysical factor is strongly dependent on the WIMP mass and annihilation cross-section in 
the latter case. 

PACS numbers: 95.35.+d, 98.35.Gi 



I. INTRODUCTION 

Flat rotation curves of spiral galaxies Q can be ex- 
plained by the presence of a dark matter halo which ex- 
tends much farther than the luminous disc. While at 
large distances the gravitational potential is completely 
dominated by the dark halo, there is still a vivacious de- 
bate about whether the dark matter is prevailing in the 
central parts of bright galaxies, and about whether its 
radial matter distribution is cuspy, or not 2\. Further- 
more, very little is known with certainty about the shape 
of the halo of disc galaxies and in many cases the halo is 
simply assumed to be spherical. 

However, if the disc is indeed an important compo- 
nent in the central parts, it should, due to its gravity, 
introduce some flattening of the dark matter distribu- 
tion. Furthermore, large scale cosmological A-body sim- 
ulations have shown that, at least at large distances from 
the center, the natural shape of dark halos is triaxial (see 
references in 0) with density axial ratios in the range 0.5 
- 0.8 (3). It is thus natural to ask whether the flatness 
and the various complex structures and substructures of 
the luminous part of the galaxy will affect the dark mat- 
ter halo, and to what extent the triaxiality of the halo 
will change a possible dark matter annihilation signal 
from the central parts of a galaxy. 

Usually, two types of asphericity are considered : flat- 
tening of the halo and departure from axisymmetry. The 
halo flattening is quantified by the value of q = c/a, 
where a is the major axis in the galactic plane and c is the 
axis perpendicular to that plane. Various observational 
methods have been used to probe the halo flattening in 
our own Galaxy and in neighboring ones (e.g. HQ). It 
was found that the measured flattening can vary over a 
wide range of values, depending on the galaxy and on 
the method used. A cross-check of the different methods 
with their systematic biases on the same galaxy would be 
welcome but is usually not possible. Measurements based 
on atomic hydrogen favour oblate halos with shortcst-to- 
longest ratios in the very wide range of 0.2 to 0.8 |fj. For 
our own galaxy and based on the thickness of the Milky 
Way's gas layer, Q argue for a rather round halo with 



flattening 0.8 < q < 1, but their result depends heavily 
on the values of the galactocentric radius R and of the 
galactic rotation speed vq 0- Recent studies of the dy- 
namics of the stellar stream coming from the disruption 
of the Sagittarius dwarf galaxy also give a wide range of 
values, between 0.5 and 1.7 [8|. 

The second type of asphericity is a departure from ax- 
isymmetry in the galactic plane. This is statistically 
quite common as a large fraction (more than 70%) of 
present day disc galaxies have bars or ovals Also, it 
is now well agreed that our Galaxy is barred in its central 
parts jl0| . Bars form naturally also in ./V-body simula- 
tions, as witnessed already in the early seventies [TTj| . 
More recently, it was realized that the presence of a dark 
halo can play an active role in the formation of the disc 
bar, if it is non rigid, i.e. if it can interact with the disc. 
Indeed, bars evolve and grow stronger by the redistribu- 
tion of angular momentum within their galaxy. This is 
emitted by near-resonant material at the inner disc and 
absorbed by near-resonant material in the outer disc and 
in the halo As a result, the halo also is deformed 

and acquires a bar structure, which is fatter and shorter 
than the disc bar, but can concern a considerable amount 
of mass [13| . 

On the observational side, departures from axisymme- 
try can be checked from the orbits of the baryons, in 
particular the HI gas that has low velocity dispersion. 
Obtaining a quantitative estimate of such asymmetries 
is, however, not trivial, since it implies a decoupling of 
the halo contribution from that of the luminous matter 
as well as a knowledge of the inclination of the galactic 
disc 01- Of course, a direct probe of the halo would 
enable to see whether the halo deformation follows the 
barred structure of the disk or not. If the recent EGRET 
diffuse gamma ray signal above 1 GeV is interpreted as 
originating from dark matter, it indeed leads to such a 
structure with an ellipticity value 0.65 ± 0.15 (l5| . 

The purpose of this article is to study the impact of an 
elliptical deformation of the halo on the expected weakly 
interacting massive particle (WIMP) dark matter anni- 
hilation signal from the galactic center. However, as 
discussed above, we do not know for certain what the 
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dark matter radial profile is, so that distributions with 
or without a cusp, and with or without a spike have to be 
considered. The density enhancement in a cuspy profile 
follows the deepening of the central gravitational poten- 
tial due to the baryon cooling through radiative processes 
[Tfil ITsl ITflj l . The presence of a super-massive black 
hole (SBH) at the galactic center can further create a 
spike, or an enhancement of the cusp, at very short dis- 
tances from the galactic center [2(J,|21| but scatterings on 
stars and capture of dark matter particle s by the SBH 
could decrease the density in this region [22, H3| ■ It is 
clear that the presence of a spike would boost the anni- 
hilation signal by several orders of magnitude. In the se- 
quel, we will see how elliptical distortions interplay with 
density profiles in the dark matter annihilation signal. 



II. HALO PARAMETERIZATION AND FLUX 
CALCULATION 

Since the main observable annihilation signal from the 
galactic center is with 7-rays |24| . we will restrict our- 
selves to this case. The observed gamma-ray flux of 
energy E, from the annihilation of dark matter parti- 
cles x (with mass m x and density p) and annihilation 
cross section UiV (into final state i), can be expressed as 
(e.g. Ref. |H) 
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where dN^/dE-y is the differential gamma spectrum per 
annihilation coming from the decay of annihilation prod- 
ucts of final state i and the integral is taken along the line 
of sight. It is customary (see (23), in order to separate 
the factors depending on astrophysics from those depend- 
ing only on particle physics, to introduce the quantity 
J(f2) for the line of sight corresponding to the di- 
rection Q, 
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We then define the astrophysical factor J(Afi) as the 
average of J(tt) over a spherical region of solid angle Af£, 
centered on the direction of the galactic center 



Iso 
Kra 
NFW 
Moore 



(3 7 a (kpc) j(l(T 3 sr) J (lCT 5 sr) 



2.0 2.0 

2.0 3.0 0.4 

1.0 3.0 1.0 

1.5 3.0 1.5 



3.5 
10.0 

20 
28.0 



2.46 x 10 1 2.47 x 10 1 

1.932 x 10 1 2.37 x 10 1 

1.21 x 10 3 1.26 x 10 4 

1.60 x 10 5 1.24 x 10 7 



TABLE I: Parameters of some widely used non spiky den- 
sity profiles models and corresponding value of J(10 _3 sr) and 
J(10- 5 sr). 
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FIG. 1: Variation of J(10~ 3 sr) with (av) for NFW and Moore 
profiles, with and without spike. The WMAP constraint on 
relic density suggests few orders of magnitude around the ver- 
tical grey line. 
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To compute the quantity J, we assume the following 
effective parameterization for the dark matter halo 



P(r) = Pa 
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A<2 



The solid angles Af2 = 10~ 3 sr and Af2 = 10~ 5 sr cor- 
respond to the angular resolutions in the EGRET exper- 
iment and the HESS and GLAST ones, respectively. The 
expected gamma-ray flux in an experiment with thresh- 
old energy Et is finally expressed as 

%{E T ) = 1-87 x KT^cm-V 1 J(Afi)Afl 



as suggested by many ACDM simulations. Rq is the 
Sun's distance to the galactic center, po is the solar neigh- 
borhood halo density and a is a characteristic length. 
The exponents a, (3 and 7 can be thought of as power 
law indices characteristic for r ~ a, r 3> a and r <C a, 
respectively. Table [I] gives the parameters for common 
halo models like the isothermal one which is not cuspy at 
the center, behaving as p(r) oc est, and the halo models 
from ACDM simulations by Kravtsov et al j^, Navarro, 
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Frenk and White (NFW) |28| and Moore et al |2g, which 
behave respectively as p(r) oc r~ 0A , p(r) oc r _1 and 
p(r) oc r at small r. Adiabatic accretion of dark mat- 
ter on the SBH could further add a central spike to these 
profiles [20j, or, if more realistic physics as off-centered 
formation of the SBH is taken into account, the SBH 
will simply enhance the cusp |2l|. Finally, scattering of 
dark matter particles by stars would substantially de- 
crease the density in the center-most region |2j| • In 
the next section we will consider both the NFW and the 
Moore profiles with and without spike. The spike's char- 
acteristic size and slope are parameterized by r sp and "f sp 
in Eq. (|SJ). We take (see 20}) r sp — 0.35 pc, and the two 
limiting values for the slope, 7 sp = 2.25 and j sp — 2.5, as 
7sp = (9 — 27) /(4 — 7). It is worth emphasizing that the 
gamma flux is dominated by the contribution from the 
unknown innermost region, especially for profiles with a 
steep slope near the center. 
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A. Annihilation effect on dark matter density 

The astrophysical factor J for the annihilation flux be- 
comes formally divergent for 7 > 1.5 (a fortiori with a 
spike) . The presence of a SBH at the center of the galaxy 
solves the problem in principle since no signal will escape 
from the region inside the Schwarzschild radius Rs (Rs ~ 
3 • 1CT 10 kpc for M BH = 2.6 • 10 6 Af Q ) . The sphere of in- 
fluence of the SBH is actually larger, as captured particles 
in a region of a few Rs are not balanced by particles scat- 
tered in from the outer shells. Following Refs.|22 | .p3 | . we 
cut the density below r c — 10~ 9 kpc: 



p(r < r c ) = 



(6) 



For halo profiles with a strong cusp/spike behavior near 
the center, the density becomes so high that the influ- 
ence of annihilations on the central density has to be 
taken into account. A simple bound can be obtained 
by letting annihilations operate in a static halo (or with 
an isotropic velocity distribution |30l|) initially with an 
infinite density, and during a time as long as the SBH 
formation time. By solving the equation 



dn / \ 2 



one obtains the upper bound 



Pn 



(av)t BH 



(7) 



(8) 



In a sense, astrophysics and particle physics aspects can 
no longer be decoupled any more in the flux calculation, 
especially for profiles with high 7 values (> 1.5). 

In Fig. ^ we can see that the presence of a spike 
near the SBH strongly enhances the annihilation sig- 
nal. The precise value of J is very sensitive to r sp , j sp 
and (o~v) /rn x . For very small values of the annihilation 
cross-section, the cut due to the capture by the SBH 
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FIG. 2: J(10 -3 sr) as a function of a) the axis ratio a/b, for 
4> ~ (blue solid line), <f> = 0.35 (= 20 deg, green long-dashed 
line) and cj> = ir/2 (red short-dashed line); b) the angle (j>, for 
a/b — 2 (blue solid line) and a/b = 9 (red dashed line). 



becomes apparent. For large values of the annihilation 
cross-section, the spike gets washed out and we recover 
the value for the profile without spike. Additional effects 
such as scattering of dark matter particles by stars or 
with baryons have been considered in the literature (|3(J), 
they compete with the annihilations for small values of 
(crv), and need to be included in a precise evaluation of 
J. 

Our main purpose, however, is to emphasize that, if 
spike or strong cusp behaviors are absent, other geomet- 
rical effects such as the asphericity of the halo become 
crucial in the J calculation, as will be shown hereafter. 
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FIG. 3: J(10~ 3 sr) as a function of prolate (blue solid line) 
and oblate (red dashed line) deformations. 
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B. Halo asphericity 

A general triaxial halo is modeled by taking 

r=[{x/af + {y/bf + {z/cf} 1/ \ (9) 

with the semi-major axis in the galactic plane aligned 
with the x axis, a > b, and abc = 1 so that the over- 
all halo mass is fixed. To study the influence of the non 
axisymmetry of the halo in the galactic plane, we set 
c = 1 and derive the variation of J as a function of a/6 
and/or the angle 4> between the x axis and the direction 
of the Sun. Note that an elliptical deformation of the 
halo can impact substantially the dark matter annihila- 
tion flux without jeopardizing the rotation curve fit. As 
shown in numerical simulations 12], the elliptical defor- 
mation is powered by the angular momentum exchange 
between the galactic bar and the dark halo. When this 
mechanism is acting solely, the deformation is strongest 
near the galactic center and decreases outwards down 
to a spherical symmetry. In principle, one could use any 
radius-dependent axis ratio function (a/b)(r) and recover 
axisymmetry at large distances (i.e. away from the bar). 
For simplicity, however, we will consider a constant ellip- 
ticity factor. This will not introduce a strong deviation 
from a realistic case because the annihilation signal is 
dominated by the contribution near the galactic center. 
Axial ratios factors of up to 5 have been obtained for the 
bar in the baryonic disk component both in observations 
and in numerical simulations, but the oval deformation of 
the dark halo is expected to be milder [Tj| . Given the un- 
certainties, however, and the unknowns from the history 
of bar formation that could well vary from one galaxy 
to another, we examine here a wider range of possible 
values. 

In Fig. |21 the variations of J(Af2 = 10~ 3 sr) as a func- 
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FIG. 4: Effects on gamma fluxes of NFW and Isothermal halo 
elliptical deformations : a) in the disc plane b) prolate and 
oblate cases . 



tion of a/b for = 0, 0.35 (= 20 deg; see review [3l| and 
references therein) and 7r/2, and as a function of 4> for 
a/b = 2 and 9 are given for the isothermal, Kravtsov, 
NFW and Moore halo profiles. For <j> — 0, J increases 
with a/b, which is expected since the higher density re- 
gion is stretched along the line of sight in this case. For 
larger values of <p, the variation becomes negative, as the 
stretching in the x axis is misaligned with the line of 
sight. Therefore, J is a decreasing function of <j> for a 
fixed value of a/b. 

It is important to notice that the impact of the ellip- 
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ticity is stronger for less cuspy halo profiles (which are 
favored by observations). Indeed, the relative contribu- 
tion to J coming from inner regions inside a small radius 
r increases with 7. Therefore, for higher values of 7, a 
larger fraction of the volume integral J p 2 dV is not af- 
fected by a change in ellipticity, as the observation solid 
angle is taken constant. 

To study the influence of the prolate-oblate shape of 
the halo, we set a = b in Eq. © and let c/a vary. The 
result is given in Fig. |31 for the isothermal, Kravtsov, 
NFW and Moore halo profiles. The variation of J with 
the prolate-oblate shape of the halo is again stronger for 
smaller values of 7. As we can see on Fig. |3 an oblate 
deformation induces an enhancement of dark matter den- 
sity along the line of sight increasing the signal, whereas 
the prolate shape decreases it and thus could be under- 
stood with arguments similar to those above. 

Finally, let's consider the most popular dark matter 
candidate, i.e. the neutralino (x), which comes from the 
neutral gauge and Higgs boson superpartners in the Min- 
imal Super symmetric Standard Model (MSSM) frame- 
work. We show in Fig. 01 the neutralino dark mat- 
ter resulting fluxes for a wide sample of supersymmetric 
models, i.e we take parameters of the MSSM to get bino 
as well as mixed bino-wino and bino-higgsino neutralino 
which have higher couplings and cross sections (see e.g. 
pf). All the points shown satisfy the WMAP require- 
ment on relic density and accelerator constraints. As the 
GLAST experiment sensitivity will probe a wide range 



of halo profiles, we clearly see that in addition to the 
(essentially inner) power law behavior of the halo, the 
geometry also alters the estimation of the fluxes and has 
to be included in flux calculations. 



III. CONCLUSIONS 

The dark matter annihilation signal from the galactic 
center has been calculated for different halo characteris- 
tics. In particular, we have shown some possible effects 
of the halo asphericity. The induced corrections are more 
relevant for flat than for cuspy cores. 

Although a plausible elliptical deformation of the dark 
matter halo does not change the expected annihilation 
signal by orders of magnitude, a consistent prediction of 
the flux from the halo shape or conversely of the halo 
shape from the signal should take those effects into ac- 
count. 
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